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We study the radiative B — >■ Kr^^'^'y decays, which are important to investigate CP violation, 
and are also relevant to assess the role of the exclusive modes induced by the & — >■ S7 transition 
to saturate the inclusive B — >■ Xs7 decay rate. Moreover, these channels do not display the same 
hierarchy as B ^ Kri^'^ modes, for which the decay into 77' is enhanced with respect to one into 
r]. The three-body radiative decays reverse the role: we find that this experimentally observed 
behaviour (although aff'ected by a large uncertainty in the case of the r;') is reproduced in the 
theoretical analysis. We compute a, B* K form factor, needed for this study, using light-cone 
QCD sum rules, and discuss a relation expected to hold in the large energy limit for the light 
meson. Finally, we examine B — >■ Rrj-y in two extensions of the standard model with universal extra 
dimensions, to obtain bounds for the compactification radii of such new physics models. 

^ PACS numbers: 12.38.Lg, 12.60.-i, 13.20. He 

CD 

Ph I. INTRODUCTION 

o 

The decay processes driven by the flavor changing neutral current (FCNC) b s transitions provide efficient tests of 
the standard model (SM) and can display deviations signaling new physics (NP) phenomena. Among such processes, 
(— I the 6 — >■ s induced decays of B mesons are the best studied and experimentally investigated; several of them have 
O I been observed through dedicated experimental analyses which have produced measurements of a variety of observables 
JL useful, on the basis of information on both inclusive and exclusive channels, to confirm the SM and constrain NP 
scenarios [T]. 

The radiative & — > 57 transition, on which we focus here, is particularly relevant. Branching fractions have been 
'— ' measured for the inclusive B — )■ Xgj mode, as well as for several exclusive channels, namely B — K*{892)^, B — 
^ i^i (1270)7, B if^ (1430)7, B K-q-f, B Kri'j, B Kcfy-f, B X*(892)7r7 and B Kinr-f, both for neutral 
^ and charged B mesons [2]. The observed exclusive modes do not saturate the inclusive rate, therefore the scrutiny of 
the exclusive transitions is mandatory in view of understanding the hadronization process for this class of channels. 
00 This is one of the motivations of the present analysis of the three-body B — >■ Krj'^'^^ modes. Moreover, there are 
other features making the multibody decays induced by 6 — 57 interesting to be studied. First, the time-dependent 
(N CP asymmetry in the neutral modes 5° — Kg j^rj^'^"/ is sensitive to NP, which may also manifest itself in producing 
right-handed photons; indeed, in the SM the photons produced in the 6 — >■ 57 transition are mainly left-handed, the 
amplitude for emitting right-handed photons being suppressed by the quark mass ratio ms/rrih [3j. Furthermore, 
the branching fractions of i? — > Krj^ and B — > K-q'^ do not obey the same hierarchy as in the two-body decays 
B — >■ Krj and B — ^ Krj\ the last process being enhanced with respect to the former one. The enhancement of 
^ two body hadronic transitions with 77' in the final state is common to several B and D decays, and is not yet fully 
understood. In the case of Dg — >• ry^'V , the gluon content of the r]' has been indicated as playing an important 

role U. For B — >■ K-q and B — )• K-q' , a possible explanation of the hierarchy between the two decay rates has been 
^ found in the destructive interference among the penguins contributions [S], and, modulo large uncertainties, this 
has been numerically reproduced in the framework of QCD factorization [5]. On the contrary, the radiative modes 
B — )■ Krj^ and B — )■ Krj'^ show the opposite trend, as one can infer from the results provided by Belle [71 18] and BaBar 
Collaborations [Sl[in], and collected in Table |l} such an outcome deserves investigations. From the experimental side. 
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Table I: Experimental results for the B — >■ Kq^'"^^ branching fractions (xlO^) from Belle and BaBar. The upper limits are at 
90% CL. 





Mode 


Belle CoUab. 




BaBar CoUab. 


B+ 


— > q-y 


8.4±1.5±JJ 


m 


7.7 ±1.0 ±0.4 tlOj 


B" 


K"qj 


8.7±i:H^I 


m 


7.1±i;J±0.4 [To] 


B+ 


^ K+q'-f 


3.6 ± 1.2 ±0.4 




1.9 ±1:1 ±0.1 (< 4.2) [9] 


B° 


^ K^q'-i 


< 6.4 


8 


l.l±i:g±0.1 (<6.6) [9] 
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the BaBar Collaboration has also measured the mixing induced (S) and direct (C) CP asymmetries in the — )■ K^rj^ 
transition. At present they are both compatible with zero: S = -0.18 ±j] j| ±0.12 and C = -0.32 ±[5:|^ ±0.07 \10\. 

The processes B — >■ Kr]^'^^ have been studied in Ref.[TT] considering exclusively the regions of the phase space where 
one of the two pseudoscalar mesons in the final state is soft, while the photon is hard. Describing the amplitudes as 
taking contributions only from virtual intermediate B* and B*, the Heavy Quark Effective Theory together with the 
light meson chiral perturbation theory (xHQET) has been employed to describe the decays in corners of the Dalitz 
plot; moreover, the 77-77' mixing has been described in the octet-singlet mixing scheme. As a result, a fraction of 
about 10% of the measured B — > A'r/7 branching ratio has been obtained. 

In the present study we improve the analysis in many respects. We take into account several possible underlying 
transitions, depicted in Fig. [l] observing that, in addition to 6 — > 57, the transition b — > sqq can contribute to the 
processes. In particular, the explicit calculation shows that the diagram (1) with intermediate virtual K* is important, 
together with diagrams (3) and (4), while the one with intermediate Ar2(1430) and the other diagrams in Fig. [I] are 
smaller. Furthermore, we do not confine ourselves to portions of the phase space, but we extend the study to the 
full three-body Dalitz plot. This requires new information, namely the B* — > K and B* ^ t] form factors, which we 
compute by light-cone QCD sum rules for a physical value of the beauty quark mass and for a wide range of four- 
momentum transferred. Finally, we consider the 77 — 77' system in the fiavor basis, in which the mixing is described 
by a single mixing angle experimentally determined with high accuracy by the KLOE Collaboration from radiative 
(f> — >■ r7^'^7 decay data. In this way, results for the SM can be obtained, together with bounds to NP extensions such 
as scenarios with universal extra dimensions (UED). 

All these topics are described in the forthcoming sections. In particular, in Sec|ll] we set up the stage of our 
calculation, considering the diagrams taken into account in the theoretical description of i? — >■ Krj'^'^^. We provide 
the expressions of the various amplitudes and identify the quantities necessary for their evaluation. Sec |III| is devoted 
to the light-cone QCD sum rule determination of the form factor which enters in the analysis oi B Krj^'^j 

transitions; we collect in the Appendix the definitions of the needed kaon light-cone distribution amplitudes (LCDA). 
The computation of the decay rates is carried out in Sec |IV[ first in the SM and then in two NP models with one and 
two universal extra dimensions. Sec fVl contains our conclusions. 

II. THE DECAYS B Krj^''>'y 

We consider the transitions B'^{p) — > K^{pk) 7;^'^(p^(/)) 7(9, e), where p, pk and are the four momenta of B, K 
and 77^'\ respectively, while q and e are the photon four momentum and polarization vector. Although we refer to the 
decays of the neutral B^ , in the following we omit the charge adopting a simpler notation; at the end of our study we 
shall comment on the charged B meson decays. The three-body transitions can be described as proceeding through 
intermediate states: The ones that we take into account are displayed in Fig. [l] The first two diagrams (1) and (2) 
take contribution from intermediate A:*(892) and K^{IA2,Q), respectively, which have width: r(A:*''(892)) = 48.7±0.8 
MeV and r(Ar2(1430)) — 98.5 ±2.9 MeV Higher kaon excitations are expected to give a smaller contribution, due 
to their larger widths and to the suppression provided by their propagators in the corresponding diagrams. The two 
diagrams (3) and (4) have B* or B* as intermediate states, which are very narrow so that we neglect their widths. The 
diagram (5) takes contribution from the intermediate 0(1020) decaying to 7^07 and having r(<?!>(1020)) = 4.26 ± 0.04 
MeV; B* is the intermediate state also in the last diagram, which involves the radiative B*B^ vertex. 




Figure 1: Some diagrams contributing to the decays B'^ — >■ K'^rj'^'^-y. The dots indicate em and strong couplings, the squares 
weak vertices. 
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To calculate the amplitudes corresponding to the diagrams in Fig. IT] we need the effective weak Hamiltonian 
describing the 6 — )■ 57 and 6 — > s gluon transition. In the SM this reads [12] : 

H,ff^~^ VtbV,: a + Cr^ O77 + CsG OsG^ • (1) 

Gj^ is the Fermi constant and Vij are elements of the Cabibbo-Kobayashi-Maskawa mixing matrix; we neglect terms 
proportional to VubV*^ since the ratio " "J^ 
written in terms of quark and gluon fields: 



is of 0(10 ). Ci are Wilson coefficients, while 0,j are local operators 



01 = {sLal'^'bLa){cLplfiCLfi) 

02 = {sLal^bLI3){CLPltiCLa) 

03 = (sLa7^&La)[(wL/37/^^t-L/i) + ••• + (^L/37M^i/3)] 

04 ^ {sLal^bLfj)[{uLplf,ULa) + ■■■ + {bLl3jf_LbLa)] 

05 = {sLal''bLa)[{uEfi7^,Uiif!) + ... + {bRfs'Jf.bRf})] (2) 
Oe = isLal'^bLfi)[{uRfiJtiURa) + ■■■ + (bRlSjfxbB^a)] 

g 

O77 = Y^[mb{sLa<y'^'''biia) +ms{sRa<y'^''bLa)]Ff^i' 



Ga 



1 zb 75 Z 

a, /3 are color indices, 5jj ^ = — - — b, and cr^'^ = 21'*''''^'^]' ^ ^^^^ electromagnetic and the strong coupling 

constant, respectively, rub and nis the beauty and the strange quark masses, while Ff^„ in Ojj and G"^, in Og^ denote 
the electromagnetic and the gluonic field strength tensors. A" are the Gell-Mann matrices. 

The Wilson coefficients appearing in ([T]) have been computed at NNLO in the standard model [I3J . The most relevant 
contribution to 6 — >■ 57 comes from the operator Oj-y, which is a magnetic penguin specific of such a transition and 
originates from the mass insertion on the external 6-quark line in the QED penguin. The term proportional to rris 
contributes much less than the one proportional to mf,, and this is the reason for which the emission of left-handed 
photons dominates over that of right handed ones in the SM. Since the coefficient Cj^ depends on the regularization 
scheme, it is convenient to consider at leading order a combination that is regularization scheme independent [14j : 

2 — 1 

where ij — ^^^-^ g^^^ C2^\^w) = 1 (the superscript (0) stays for leading log approximation); furthermore, 

14 16 6 12 

"^"23 "'"23 "^'23 "*-"23 

as = 0.4086 ag = -0.4230 07 = -0.8994 ag = 0.1456 

/ii = 2.2996 /i2 = -1-0880 ^3 = ^4 = -^ (4) 

= -0.6494 he = -0.0380 = -0.0185 hg = -0.0057 . 

The effective weak vertex Oj-y contributes to the diagrams (1-4) in Fig. [l] through hadronic matrix elements that we 
define below. However, before doing that, we turn to the rj — rj' system. 

The r] — rj' mixing is usually described in two different schemes, adopting either the singlet-octet (SO) or the quark 
flavor (QF) basis, and in each scheme two mixing angles are involved [T3]. Here we adopt the quark flavor basis 
defining 

M - -^{\uu) + \dd)) 

\Vs) = \ss} , (5) 
so that the 77-77' system can be described in terms of the mixing angles ipq and ip^: 

= cos ipq\r]g) - smifislVs) 
\r]') = sin(fq\r]q) + cosips\iis} ■ (6) 



4 



The difference between tpq and ips is due to OZI- violating effects and is experimentally found to be small {(pq — (ps < 5°), 
so that it has been proposed that the approximation of describing the rj — rj' mixing in the QF basis and a single 
mixing angle is suitable [TS]. The simplification ~ (ps ^ f is supported by a QCD sum rule analysis of the 
4> ^ r]"f and (p — > r]'"f decays |16| . A precise determination of the rj — r]' mixing angle has been obtained by the 

KLOE Collaboration measuring the ratio — ^ in the flavor basis with a single mixing angle, with the result: 

r(0 777) 

(f = (41.5 ± O.Sstat ± O.Tsysi ± 0.6th) [II]- This analysis has been improved performing a global fit of the transitions 
V — >■ P7 and P — >■ Vj {V = 4>, p and P — tt^, 77, 77'), allowing a gluonium content in the 77' and including the 

r(^' — ^ tt) 

measurement of the ratio — ^-7; r 118 . The outcome is that, even though the gluonium content of the 77' is 

r(7r"' — > 77) 

significant, the result for the 77 — 77' mixing angle is only negligibly affected. Therefore, we set ip to the value quoted 
above. 

Let us now consider in turn the various diagrams in Fig. [l] 
• Diagrams 1 and 2 
The corresponding amplitudes read: 



A2 = A{B^K*-i) 



s - mj^, + i niK' Tk' 



AiK* Kt] 



(7) 
(8) 



with 



A{B Kl^^-f) = Ce*'^ ^(rTifc + 771,) < Kl2){pK,l)\s(T^,^q''h\B{p) > +{mb - rris) < K*^2)iPK,e)\s<7^,^q''-f5b\B{p) >J , 

to be computed for an on-shell {q^ — 0) photon, defining s ~ {p — qY = Af^^^,,. The factor C is C = 

A—j=VthVi^c'f^^^ ^^ ^ . The hadronic matrix elements of weak Hamiltonian operators are parametrized in terms 
of form factors: 



<Kl2){pK,i)\-s<7^Xh\B{p)> = 7e^,„0?* 2 Ti '-'■"{q- 
< K(2){PK,e)\s(T^^,q''j5b\B{p) > = ~e*^{Ml - s) - {e* ■ q)(p + pk) 



9m 



Mi 



(9) 



(10) 



with e denoting the polarization vector of the K^^^^^ mesons; in the case of iir2(1430), which is a spin 2 particle, the 
polarization is described by a two indices symmetric and traceless tensor, therefore in (9p0 1 it is understood that 



-gg PP 



The condition holds: T, 



(0) = T2 



' (0) . The variable s in the definition of the hadronic matrix 



elements takes into account that the ™6sons are off-shell, and is needed to ensure gauge invariant amplitudes. 
In the same diagrams strong vertices also appear, which are defined as follows: 



A{K*^Krj^>) = gK'Kr,(')^-P^in 



(11) 

(12) 



Within the fiavor scheme for the 77—77' mixing the relations gx^Kr] = (cos ip+y/2s'm p)gK'+K+Tr" ^-nd gx-'Kr]' — {s'mip — 
\/2 cos 'f)gK*+K+TT° can be worked out. Assuming the width of K*^ saturated by the two modes K*^ — > /f+7r°, K'^tt^ , 
and using the relation gK*+K°'K+ = ^/^gK*+K+-K'>^ from r(A'*+) = 50.8 ± 0.9 MeV we obtain gK*+K+-K^ = 6.5 ± 0.06. 

The strong coupling gK*Kri can be estimated, although with a large uncertainty, using the measurement 'S>{K2 — > 
Krf) = (1.5±f:;^) X 10-3 P] together with r(Ar;^), obtaining: gx^Kn = 1.43 ± 1.60 GeV-^ On the other hand, no 
information is available for gx^Kri'] however, since, as we shall see, the contribution of this diagram is small in the 
case of 77, it is reasonable to neglect it also in the case of the 77' in the final state. 
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• Diagrams 3 and 4 
The two amplitudes read: 



with 



and 



^3 = A{B ^ B*rj'^''>) 
Ai = A{B B*K 



^A{B* ^ Kj) 

t — rrig, 

('). 



u — nif 



-A{B: ^ ,7^'^ 7) , 



(13) 
(14) 



AiB*^K-f) = Ce*''[{mb + m,)<K{pK)\sa^,q-'b\B*{p\e)>+{mb-m,)<K{pK)\sa^,q''^5b\B*{p',e)>] , 



A(b: ^ r,W7) = Ce*" (mj, + to,) < r;('^(p^(„)|sa^.q''6|B;(p',e~) > +(mb - m,) < r^^'^ip^u,)\sa^,q''j5b\B:{p\i) > 



(')( 



<K{pK)\-sa^,q^h\B*(p',i)> = *e.ra/^?V>^ 2 Tf '^^(g^) 
< K{pK)\sa^,^q^'-izh\B*{p\~e)> = U^{t - m\) - {i ■ q){p' + px),. 



t — m\ 



{p' +Pk)i. 



(15) 



(16) 



<V^'Hp,^i,,)\sa,,,q^'b\B:ip',~e)> = ie^^^p^p'^p" 2T^ 



(') 



< V^'HPrj('))\s(T^^q''j5b\Bl{p',e) > ei.(w - m^co ) - (? ■ g)(p' + P^w )^ 



Tf--''\q^) 



1— (P' +^',,('))^ 



rjV) 



(9^ 



(17) 



(18) 



(p'j e) denote the four momentum and the polarization vector of the -B^*^^; moreover ^ (0) — T2 ^ (0) and the 

same for ri^'\ The variable t — (q + p^)^ takes into account the ofF-shellness of the B* in diagram (3), while the 
variable u = (p^(/) +q)^ accounts for the off-shellness of the B* in diagram (4); obviously, s + t + u = M'^ + m\ + rn^^iy 
As for the strong vertices appearing in the two amplitudes, we define 



A{B B*r]^'^) 
A{B B*K) 



9B'B'n<o e -Prji'-) 
Qb'BK ^* ■ Pk ■ 



(19) 
(20) 



The two couplings 5b->b,,(') and gs^BK can be obtained, invoking SU(3)f symmetry, from the analogous quantity 
gB-B-TT- 9B'Bri = costpgB'BTj, = cos(pgB'ByT, gB'B-qi = sin'^gB-Br^^ = sim^gs-B^ and gB^BK = gB'B^- As for 
gB-'BTTi it can be related to a low-energy parameter g that describes the coupling of heavy mesons belonging to the 
doublet of heavy-light quark states with spin-parity = (0^, 1^) to light pseudoscalar states in the framework of 
the xHQET [19]: gB*Bn = ^^^5- There are several theoretical determinations of g spanning the range [0.2, 0.5] [20] . 
However, g can be extracted from the measured decay width of D*^ — D^tt"'" [21], obtaining g — 0.59 ± 0.01 ± 0.07 
[22]. We use this value in our analysis. 

• Diagram 5 

The contribution of the intermediate 0(1020) is represented by the amplitude 



A5 = A{B Kcj)) 



Adopting factorization, the first amplitude in (21) can be written as 

A{B -> K<P) = ^VtbVtMKipK)\sjf.{l - l5)b\Bip)) (0(p^, 6~)|s7^s|0) 



(21) 



(22) 
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where a^, is an effective Wilson coefficient that we set to the value = 0.064 ±0.009 from the experimental branching 
fraction 'B(B ^ K (jj) = {S.6±\-l) x lO'^ 0. Furthermore, we use the parametrizations 

{K{pK)\-S7^{l~l5)b\B{p)) = .f^^''{q'){p + PK), + f^^''{q'){p~PK)^ 

e)|s7ps|0) = /^m^?* . (23) 



When these two definitions are inserted in (22) only the form factor ff^^ contributes, and we adopt for it the 
determination in Ref.[53]. The value = (0.232 ± 0.002) GeV comes from the experimental datum 'B((/) — >■ e~^e~) = 
(2.954 ± 0.030) X 10"''. 

Following Ref.[T5], the amplitudes yl((/) rj^'^^) can be written as 

A{4> ^ ri^'h) = -lF^^^"'^{q^)e,^pse*^{p^r{p,,)f'e' . (24) 

The form factors i^"*"^'' ^(q^) were determined using QCD sum rules, providing their values at = (multiplied 

by the strange quark charge in units of e): \g*'i^\ = -F't'^'^^{Q) = (0.66 ± 0.06) GeV-^ and = ^^'^^'''''(O) = 
3 3 

(1.0 ± 0.2) GeV""'^ [1^, results used in our analysis. 
• Diagram 6 

It is possible to show that the diagram (6) provides a tiny contribution with respect to the others. 
Let us discuss this in the case of the rj. The amplitude can be written in terms of A{B* — > B^) and 
A{B* — Kvi). In order to understand how large this contribution is, we can invoke naive factorization, writ- 
ing A{B* ^ Kri) = ^V*^VusaJ^ {K\s^i'{l - -iz)\B*) {■q\u'y ^{l - 75)u|0), with aj^ ~ -0.286 an effective Wilson 
coefhcient for color suppressed decays. The ry-current-vacuum matrix element involves (in the flavor basis for the 
•q — rj' mixing) the constant — fqCOS(f> with fq ~ (77|m7^(1 — 75)m|0) = ^ (Pj;)^i- On the other hand, the 
matrix element (i^|s7^(l — 75)|i3*) can be decomposed in terms of several form factors; however, when contracted 
with {pri)^n only one of such form factors contributes, usually denoted as ^o(w^), which in the large energy limit 
of the final light meson coincides with computed in the next section. The other ingredient is the radiative 



amplitude A{B* B^), which can be written as A{B*{p' , e) B{p)j{q, e)) = e - — I- -p- ^afira e*" e^p'^p'"" , with 

e-b (sq) the b (q — d) quark charge in units of e. A determination of the mass parameters A;, and A^ can be found in 
|24j : Ab = 4.93 GeV (close to the b quark mass) and A^ — 0.59 GeV. As a result, the contribution of the diagram (6) 
to the branching fraction is O(10~^'^). Therefore, in the following we neglect this amplitude. 

As it emerges from the above discussion, important quantities are the form factors appearing in the diagrams 
(l)-(4). In the next section we compute ~^^{<l^) by hght-cone QCD sum rules [5S]. SU(3)f symmetry and 

the QF 7? - ?]' mixing scheme allow also to fix: Tf '^''(0) = -sunpTf^^ (Q) and rf=""''(0) = cos Tf *^^(0). 
As for Ti^^ , several determinations can be found in the literature; we use the light-cone QCD sum rule result 
^^*(0) = 0.333 ± 0.028 [21, to be consistent with the determination of Tf'^^. This value is compatible with 

the one obtained by three-point QCD sum rules based on the short-distance expansion [271. Finally, for ^ we 

use rf (0) = 0.17 ± 0.03 ± 0.04 [2g. 

There is a remark concerning the relative strong phases among the various amplitudes. While the sign between 
the amplitudes (3) and (4) can be fixed invoking ^HQET and the flavour symmetry, the relative phase between, e.g., 
(1) and (3) does not follow from theoretical arguments. Therefore, we consider it as a parameter to be determined 
empirically from the experimental data. The phases appearing in the other amplitudes do not play a role in the 
branching ratio due to the small size of such diagrams. 



III. FORM FACTOR Tf ~*''((7^) BY LIGHT-CONE QCD SUM RULES 

To compute the form factor {(f') by light-cone QCD sum rules (LCSR) we consider the two-point correlation 

function with the external kaon state 

U^.u{p'.q)=^ J d'xe"^-^K{p')\T{J^{x),K{0)}\0) (25) 
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where = scja^q^h is the quark current appearing in the matrix element (15). = b^^q is the vector current with 
the quantum numbers of the B* meson {q = u, d), and its matrix element between the vacuum and the B* state is 
parametrized in terms of the decay constant Jb* , 



{B*{p' + qre)\h.q\Q) = fB-mB'K 



(26) 



The LCSR method consists in expressing the correlation function Eq.(25l both in QCD and in terms of a hadronic 
representation. can be decomposed in independent Lorentz structures, one of which can be used to compute 



n,ti^(p', q) = * ^tiaTiyq^p'^^iiv' + qf-, q^) + other structures 



(27) 



In terms of hadronic states, the correlation function in (25) can be written as 

nHAD. , , {K{p')\J^\B*{p' + qre)){B*{p' + qre)\VM , sr-{K{j>')\J^\h{p' +q)){h{p' + q)\V,\Q) 



to|. - {p' + qf 



and consists in the contribution of the B* meson and of the higher resonances and of the continuum of states h. The 
first term in (28) contributes to the invariant function ll{{p' + q)'^,q'^), since 



{Kip')\J^\B*{p' + q,i)){B*{p' + q, e)|K|0) 



I e 



2Tfy^iq^)fB.mB> 



(28) 



m^, — [p' + q)'^ nig, — (p' + qY 

In a one-resonance+continuum formulation, the hadronic representation of the function Ii{{p' + q)^ , q^) can be written 



as 



n 



HAD 



{{p' + qf,q^) = 



2T{ 



{q^) fB*mB* 



rn%, - [p' + qf 



s - {p' + qf 



(29) 



So 



where higher resonances and the continuum of states are described in terms of the spectral function p''^{s,q^) which 
contributes starting from a threshold sq. 

The QCD expression of the correlation function is given by 



^Qc°((p' + ,)^q^ 



as 



{p' + qf 



(30) 



This expression comes from an operator product expansion (OPE) of the T-product in Eq.(25) on the light-cone, 
which produces a series of operators, ordered by increasing twist, the matrix elements of which between the vacuum 
and the K (required to evaluate Eq.(25 )) are parametrized in terms of K light-cone distribution amplitudes (LCDAs). 
The equality of the hadronic and QCD representations of the correlation function, Eqs.(29) and (30), does not yet 
allow us to derive the Ti form factor, since the hadronic spectral function p'* is unknown. However, we can 

invoke global quark- hadron duality above the threshold Sq [2§] , which amounts to identify integrals of the spectral 
function p'' with corresponding integrals of pQ^^ — ^Imll'^'-^^ , and in particular 



ds- 



p\s,q^) 
-{p' + qf 



ds 



imnQCD(5^^2) 

s ~ {p' + qf 



(31) 



Using global duality, together with the equality n^^°(p',g) — H'^^^ {p' , q) , from Eqs.(29) and (30) the equation 
follows: 



2T/^ ^^(g^)/B»mB» 

- [p' + qf 



so ImnQCD(s^g2) 

ds - 



{p' + qf 



(32) 



The subtraction of the continuum and of the higher-twist contributions, leading to (32), can be optimized, following 
the QCD sum rule procedure, by a Borel transformation of the hadronic and of the QCD expressions of the correlation 
function, hence of the two sides in Eq.(32). This transformation, which applied to a function 3'{Q^) (with Q"^ = — g^) 
is defined as 



®[J(Q2)] ^ li 



>oo. n^oo. 



1 

(n- 1)! 



(-Q2)" 



d 



HQ') 



(33) 
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where is the Borel parameter, produces the equahty 



1 



(S + Q2)n 



cxp(-s/M2) 



(34) 



This operation improves the convergence of the OPE series by factorials of the power n, and for suitably chosen 
values of AP enhances the contribution of the low lying states to the hadronic expression of the correlation function. 
Applying the Borel transformation to both n^^^ and n'^'-'^ we obtain: 



2rr-^(,2)/^.^^. exp(^-^)=i |"dsexp(-^) lu^Ii^^^ [s , q') 



(35) 



The calculation of II'^'-^'-', based on the expansion of the T- product in (25 1 near the light-cone, involves matrix elements 
of nonlocal quark-giuon operators. The final sum rule for Ti has the form: 



2rf -^(g2) fg.ruB' = fiQCD(o) ^f^QCD(i) 



(36) 



where the symbol H indicates that Borel transformation and the continuum subtraction have been performed. n'SC'^C) 
gets contribution only from two-particle distribution amplitudes, while is written in terms of the three-particle 

ones, all collected in the Appendix. Their expressions are: 



du 



u 



mb(t)Kiu) + 



K 



IvP , 

1 



1 du 



*4/f(w) 



U(pp{u) + -(j)a{u) 

6 



{so - q") 



(37) 



— 04k(u) e 



+ 



(m^ — q^y- 



So 



AP 



"o04k("o) 



jjQCD(l) 



da I 



1 — ai 



da? 



dv e 



3K- 



{ai,l - ai - as^as) {ml-q^)\ ^ uibf, 



$(q;i, 1 - ai - as, as 



AP 

2 „2\ 



+ -jjii^vAo'i, l~ai~ as, as) 



/ (mg^\ 
' \ 2C/M2 J 



e M 

+mbjK 7:^2 — 



{mi - q2) 

Ug 

da I 



"3 
das 

1(0 — Qi ^-^3 



-l>(uo, 1 - Wo - as, as) 



da I 



-uo - ai 



-l>(ai,l - uo,uo - ai) 



(38) 



Wo - ai rub/K _ / -, X 
J3K (P3K[ai, 1 - ai - as, as) H ^ 9'P±{(^ii 1 - ai - as, as) 



Upas 



+ 



rribfk 



^(ai, 1 - ai - as, as) 



M2 

+mbfKe 



da:. 



l-un 



dv V ( 



^(«3) A 



^(l-^o) 

Uo(l - Uo) 



^vl,(as) 



l-tio "3 
2 



2wAP 

1 - Mo 
li2^/2 



In the previous equations we have defined uq = — U ^ai-t-wasjUi^l — was. Furthermore, the LCDAs have 

So - r 
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been combined as follows: 

$(ai, 1 - ai - as, as) 

1 - ai - as, as) 

^(aa) 
1 - ai - as, as) 
l>(ai, 1 - ai - as, as) 
#(as) 



[<p±(t, 1 - < - as, as) + (^11 (t, 1 - t - as, as)] 
- t,0,i) 

dt [ip±{t, 1 - < - as, as) + (^11 (t, 1 - t - as, as)] 



tpj_{ai, 1 - ai - as, as) - (fj_{ai,l - ai - as, as) 
<l>(ai, 1 - ai - as, as) - $(ai, 1 - ai - as, as) 
*(as)-^'(as) . 



(39) 



The distribution amplitudes entering in the previous relations can be classified according to the twist: (j)K is a 
distribution amplitude of twist 2, 0p, and (p^K are twist 3, (piK, V'4-R': f±j 'Ph 'f'-^ twist four. We have 
used the definitions of various matrix elements defining the LCDAs as well as the updated numerical values for their 
parameters in Ref.(3D]. We use the value Mb- — 5.325 GeV, and the quark mass mf, = 4.8 GeV which also enters in 
the calculation of the SM effective Hamiltonian in the next section. 

Eqs.(36),(38l and (39) allow us to compute once the threshold sq and the Borel parameter Af^ are fixed. 



The threshold is set to sq = (33 ± 2) GeV , a value appearing also in QCD sum rules involving the B* meson [31], 
which is close to the estimated mass squared of the first radial excitation of B* . For each value of squared momentum 
transfer , the form factor also depends of the Borel parameter , which can be fixed requiring stability 

against variations of M^. In Fig. [2] we depict the dependence of the Tf ~*'^(0) on the Borel parameter M^. The 
band reflects the uncertainties on the other quantities entering in the calculation, including the uncertainties on the 
LCDAs parameters quoted in |30| . on the threshold sq and on fs- for which we use the value computed by QCD sum 
rules in 'ST]: fs' — 0.195 ± 0.035 GeV. Although the results are quite stable with AP, in the numerical analysis we 
fix the stability window in the range = (17 ± 3) GeV^. 

0.5 



0.4 
0.3 
0.2 
0.1 
0.0 



10 12 14 16 
(GeV^) 



18 



20 



Figure 2: Dependence of Ti ~^ (0) on the Borel parameter M . 

For all values of in the range [0, 20] GeV^ the computed form factor is plotted in Fig. [3j The functional 
dependence is obtained fitting the sum rule result by a single pole parametrization: 



(0) 



(40) 



1 - 



Ml 



with 



Tf*^^(0) =0.30 ±0.066 



Mp = 5.767 GeV 



(41) 



The uncertainty of ^^(0) also accounts for the variation of the Borel parameter within the stability window. 

We conclude this section with a comment about the relations among the form factors that parametrize the B to 
a light meson L matrix elements. As shown in Ref.|32j, when the energy of L in the rest-frame of the decaying B 
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meson is large, the form factors describing the B to L transitions can be related among each other. Considering 
also the heavy quark limit, it is possible to relate the B ^ L form factors to the B* ^ L ones, as shown in [33] 
where also the perturbative corrections to the large energy relations have been worked out. In particular, the relation 
— ff~*^ should hold for large kaon energy, where /+ parametri zes the matrix element of the quark vector 
current 57^6 between the kaon and the B meson as in the first equation in ( ps] ). Since large kaon energy means q'^ close 
to q'^ = 0, one might exploit the relation ~*^^(0) = /+~*'^(0), as done in [TT]. Our computed form factor 
together with ff^^{Q) = 0.331 ± 0.041 determined by LCSR ^3], fulfills the relation within the uncertainties. 



IV. B Kri'^'^-y DECAY RATES AND PHOTON SPECTRA 

In Secs|II|and |III| we have collected the quantities necessary to evaluate the amplitudes corresponding to the diagrams 
in Fig. [1] Input parameters are the quark masses, mt already fixed and nis — 0.130 GeV [53], and the CKM matrix 
elements Vtb = 0.99, Vts = 0.04 [2]. We do not include the uncertainties on these parameters because they are small 
with respect to the uncertainties of the other input quantities; moreover, mg plays a negligible role in the final result. 
We set the renormalization scale at which we compute the coefficient c\tj^'' to /i = 5 GeV, close to the b mass. 

We can now discuss the branching fraction and the photon spectrum of i? — > if 777 and B K-q'^ in the SM and 
in two new physics models with universal extra dimensions described below. We anticipate that such new physics 
scenarios belong to the class of minimal flavor violation models, therefore the only modification with respect to the 
SM consists in a different value of the Wilson coefficients in the effective weak Hamiltonian. Therefore, the three 
cases, the SM and the two UEDs, share common features, namely the hierarchy among the various decay amplitudes 
and the shape of the photon spectrum. 

For what concerns the intermediate states in the B — > Krj'^'^j decay amplitudes, the most important contributions 
are represented by the diagrams (1), (3) and (4) in the case of the 77, while for rj' the first diagram contributes much 
less than diagrams (3) and (4). This is due to the coupling gK'Kri which is much larger than gK'Krj''- indeed, from the 
relations in Secjiljwe get gK'Ki) = 11 ± 0.1 and gx-Kr)' = — 2.57± 0.19. We only consider a phase 9 between the sum 
of the first two amplitudes and A3 + A4; the fifth diagram turns out to be much smaller than the others, hence we 
assign to it the same phase as Aj, and since a wide change of its phase does not modify the result. In the case of 
77' we do not consider the contribution of if 2 : ^-^id 9 is the phase between Ai and A3 + A4 + . From the calculation 
of !B(B — >■ Krjj) we shall see that there is a range of values of 9 allowing to reproduce experimental data in Table |l] 
Let us start from the standard model. 



A. B — > Kri'^'^'-f in the standard model 



The plot of the computed 23(5 —> Krjj) as a function of the strong phase 9 is depicted in Fig. B] The experimental 
results in Table |l] can be obtained for 9 = 1.8 ± 1.0 rad, corresponding approximately to 9 = (7 ± 3) (we have 
considered the range < 9 < tt since the plot is symmetric with respect to 9 = 0; indeed, in the branching ratio 
the term proportional to sin 9 is two orders of magnitude smaller than the one proportional to cos 9). For the central 
value of 9, the photon spectrum is depicted in Fig. [5] It is peaked at large photon energies and has a structure as the 
effect of the virtual K* . The Dalitz plot in the plane {M,jk, Eri), displayed in Fig. |6j also shows the effect of the K* 
in B — > if 777, at the limit of the phase space: it should be observed in the data. 
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Figure 4: Branching ratio 3(_B — )■ Krj'y) as a function of the strong phase phase 6 between Ai + A2 and + A4, + A^,. The 
horizontal band corresponds to the experimental result. 



In the case of the 7/', since the diagram (1) gives a small contribution with respect to (3) and (4), strong phases 
do not play any role. The result for the branching ratio is 23(5 — >■ Krj'j) = (2.78 ± 1.14) x 10~^, with the photon 
spectrum depicted in FigjS] and the Dalitz plot shown in Fig|6] The theoretical result of the branching fraction for 
the neutral mode is compatible with the experimental datum in Table |Tj which is affected by a large uncertainty. 
The experimental error is smaller in the charged mode: in this case, while the BaBar result is compatible with the 
calculation, the Belle measurement is larger. Before commenting on the charged mode, it is worth observing that, for 
the B — >■ Kr/^'^j three-body channels, the hierarchy between the modes with rj and 77', observed in data, is reproduced 
by the theoretical calculation in the frameworks of the QF scheme for the rj ~ rj' mixing. 
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Figure 5: Photon spectrum in B° — >• K^rij (left) and — ^ K^'rj'j (right). The phase 9 is set to 9 — 1.8 rad. 

Let us discuss the differences between the neutral and the charged B^ radiative decays. The analysis of the 
charged modes would be the same as the one we have presented for the neutral modes, except for the contribution 
of the inner bremsstrahlung (IB) diagrams with the photon coupled to the charged initial B~^ and final mesons. 
The kinematical region in which the bremsstrahlung could be competitive with the other decay mechanisms is for soft 
photons, due to the presence of a pole at vanishing photon energy. To describe this contribution to the 5+ K'^rj'^'^^ 
decay, we invoke the Low theorem ^34j which, for scalar particles, allows to relate the amplitude of the radiative mode 
with a soft photon to the amplitude A{B^ — > A'+t]*^'^): 



AiB {B+ {p) ^ K+ [pkW'^ (p^<„ )7((7, e)) - e 



e- -PK e -Prii') 



q-PK 



q ■ PrjC) 



A{B+ K+ri 



(0) 



(42) 



The two-body amplitudes in (42) can be determined from the experimental branching fractions !B(i?+ — > K^if) = 
(2.33±g:i)xl0-^ and 2(5+ -> K+r]') = (7.06±0.25) x 10"^ [2]. In the case of the 7;, the estimated contribution of the 
IB diagram to the decay rate is of order 0(10^^*); therefore, the rate of the charged mode is not significantly affected 
by this effect, as indeed observed in the data in Table |l] The contribution is more important for the rj': varying the 
relative phase Oi,rem between the bremsstrahlung amplitude and the other considered amplitudes, 23(5+ — > K^rj'^) 
varies between (2.3 ± 1.0) x 10"'^ and (3.8 ± 1.4) x lO"'^. 



12 




3%(GeV) M^-jr(GeV) 

Figure 6: Dalitz plot of 5° li^r)^ (left) and W K^ri'^y (right) in the plane (M^(,);f , -E,j(,)). The photon and the 77''^ 

energy in the B rest frame are i?-y = (^M% — Af^(,)^^ /2Mb and -E^(/) = (^M% + — /2Mb, respectively. Lighter colors 
correspond to higher values of the distribution. 



B. B — > Kr]-^ in the ACD new physics scenario 



Now we consider a new physics scenario: The Appelquist-Cheng-Dobrescu (ACD) model with a single universal 
extra dimension (UED) 35 . The model is a minimal extension of SM in 4 + 1 dimensions, with the extra dimension 
compactified to the orbifold IZ2 and the fifth coordinate y running from to 2t:R, y = and y — ttR being fixed 
points of the orbifold. All the fields propagate in all 4 + 1 dimensions, therefore the model belongs to the class of 
universal extra dimension scenarios; one of its motivations is that it naturally provides candidates for the dark matter. 

In the ACD model the SM particles correspond to the zero modes of fields propagating in the compactified extra 
dimension. In addition to the zero modes, towers of Kaluza-Klein (KK) excitations are predicted to exist, correspond- 
ing to the higher modes of the fields in the extra dimension; such fields are imposed to be even under a y — > —y 
transformation in the fifth coordinate. On the other hand, fields which are odd under this transformation propagate 
in the extra dimension without zero modes, and correspond to particles without SM partners. 

The masses of KK particles depend on the radius R of the compactified extra dimension, the new parameter with 
respect to SM. For example, the masses of the KK bosonic modes are given by 

ml = ml + j^ n = l,2,... (43) 

Too being the mass of the zero mode, so that for small values of R, i.e. at large compactification scales, the KK 
particles decouple from the low energy regime. Another property of the ACD model is the conservation of the KK 
parity (— l)-*, j being the KK number. KK parity conservation implies the absence of tree level contributions of 
Kaluza Klein states to processes taking place at low energy, forbidding the production of a single KK particle off 
the interaction of standard particles. This permits to use the electroweak measurements to provide a lower bound 
to the compactification scale: l/R > 250 — 300 GeV 36^. Moreover, this suggests the possibility that the lightest 
KK particles, namely the n — 1 Kaluza-Klein excitations of the photon and neutrinos, are among the dark matter 
components [571155] . 

Since KK modes affect the loop-induced processes, fiavor changing neutral current transitions can constrain this new 
physics scenario. Indeed, many observables are sensitive to the compactification radius in the case, e.g., of processes 
involving B, Bg and Af, [5^43) . In the ACD model no operators other than those in ^ contribute to 6 — >■ 37, and 
the effects beyond SM are only encoded in the Wilson coefficients of the effective Hamiltonian. The contribution of 
KK excitations modifies in particular the coefficient Cr-y, which acquires a dependence on the compactification scale 
l/R. For large values of l/R, due to decoupling of massive KK states, the coefficient Crj (whose explicit expression 
can be found in [55] ) reproduces the standard model value. 



In Fig. 7(a) we compare the theoretical prediction for 23(_B —> Kr]^), obtained in this framework, to the experimental 
datum, having fixed the phase 6 — 1.8 rad. The uncertainty on the theoretical result does not allow to constrain l/R 

in a stringent way, since only the bound — > 250 GeV can be established, lower than the bound provided by other 

R 
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processes [3D]. 



C. B — > Kr]'y in a model with two UEDs 

Models with universal extra dimensions may have a larger number of additional dimensions, and indeed a model 
with two compactified UEDs has been proposed in [44] . The two extra dimensions are flat and compactified on a (so- 
called chiral) square of side L: < x'^,x^ < L, where and are the fifth and sixth extra spatial coordinates. The 
compactifi-cation is performed identifying two pairs of adjacent sides of the square: (y,0) — (0,y) and {y,L) — {L,y), 
for all y € [0,L], which amounts to folding the square along a diagonal. The fields are decomposed in Fourier modes 
in terms of effective four dimensional fields labeled by two indices {I, k). Hence, the KK modes are identified by two 
KK numbers which determine their mass in four dimensions: zero modes corresponds to SM fields. The values of the 
fields in the points identified through the folding are related by a symmetry transformation. For example, for a scalar 
field, the field values may differ by a phase. The choice of the folding boundary conditions (and of the constraints on 
such phases) is mostly important in the case of fermions, since a suitable choice allows to obtain chiral zero modes, 

while higher KK modes have masses determined (as for scalars) by the relation: Mi^k = , where i? = — is 

the compactification radius. The theory has an additional symmetry, the invariance under refiection with respect to 
the center of the square. Such a symmetry distinguishes between the various KK excitations of a given particle. A 
KK mode identified by the pair {I, k) of indices changes sign under reflection if Z + A: is odd, while it remains invariant 
if / + fc is even. As a consequence, the stability of the lightest KK modes is guaranteed, and such modes are good 
candidates for dark matter. 

The model comprises the SM particles and their KK excitations, together with new particles without a SM corre- 
spondent, described by fields whose Fourier decomposition does not contain a zero mode. Examples are the mixing 
of the fifth and fourth components of the vector fields and the Higgs fields. All these new particles may contribute as 
intermediate states in the FCNC loop diagrams and, as in the single UED case, they modify the values of the Wilson 
coefficient in the effective Hamiltonian o| without introducing new operators. The explicit expression of C^!^^ in 
this model can be found in Ref. 45 . It should only be mentioned that the sums over the KK modes entering in the 
expression of the Wilson coefficients in the extra dimensional framework diverge logarithmically, and should be cut 
in correspondence of some values of Nkk = I + k, viewing this theory as an effective one valid up to a some higher 
scale. The condition Nkk — 10 has been chosen in [45] . 

With the modified value of C^lj^ , which depends on the compactification radius, we obtain the B — > Krj^ branching 

versus 1/R. Comparison with experiment provides the constraint — > 400 GeV, which 
■ I R 
is close to the bound established from other processes [HI US] . 



fraction plotted in Fig. 



7(b) 
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Figure 7: Branching ratio S(B — >■ Krjj) as a function of the inverse of the compactification radius (in GeV) in models with 
a single UED 7(a) and two universal extra dimensions 7(b) The strong phase 9 is set as in the SM. The horizontal band 
corresponds to the experimental measurement. 



14 



V. CONCLUSIONS 



B decays to two light pseudoscalar mesons and a photon are interesting, as witnessed by the experimental efforts to 
determine, for the modes B — > Krj^'^^j considered here, the branching fractions and the CP asymmetry parameters. We 
have studied such channels considering the contribution of amplitudes corresponding to several intermediate states, 
K*(d,Q2) and 1^1(1430), as well as B*, B* and ^(1020). A Ught-cone sum rule determination of the form factor 
Ti ^^{q^) has been performed: this form factor is of interest since it also enters in other amphtudes involving B* 
mesons. Introducing a strong phase 6 between the the first two considered contributions and the other three, we have 
shown that the measured B — >■ Krjj branching fraction can be reproduced. On the other hand, the experimental 
uncertainties in !B(i3 K'q'-y) are large, so that the comparison with the theoretical result does not provide constraints, 
at present. In any case, the modes with r]' in the final state are not enhanced with respect to those with the r], 
as experimentally observed. The photon spectrum, as well as the Dalitz plots, are sensitive to the intermediate 
contributions. 

We have studied the same transitions in NP scenarios with one and two universal extra dimensions. In particular, 

we have found that, in the case of two UEDs compactified on the chiral square, the bound — > 400 GeV can be 

H. 

established from 'B{B ^ Krj^). 
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Appendix: LCDAs of the K meson 



Here we collect the matrix element defining the LCDAs of the kaon used for the calculation of the form factor 
• Two-particle LCDAs 



Jo ^ P ' X Jq 



{K(p')\-s{x)i^^qmQ) 



mq Jo 



rus + nig 



{K{p')\s{x)a^„^5q{0)\0) = i {p'^x^ - pl^x^j) 



fKm 



K 



• Three-particle LCDAs 



Q{ms + niq) Jq 



{K{p')\-s{x)g,G'-ap,^^qm{)) = 



ihKWp'pQl-p'^p'X) - {p'^'X -P^PJp)] / D«.03/f(Qi,a2,a3)e 



{K{p'Mx)gsG'^j^j,q{0)\0) 



Ik 



J Daie*f'--("i+''"^)<^i(ai,a2,a3) 



p' ■ X J " 
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Ik 



P 9 



X p 



rip 



p' ■ X 



■P \9 



rip 



+L 



K 



P' 



rip 



{p'^X^ 



p' ■ X 



p X 



p' ■ X 



The following definitions have been used: T)ai = daida2da3 5(1— ai — a2 — as) and G"^ = \e°'^^'^GeT- The expressions 
of the LCDAs listed above, together with the numerical values of the parameters entering in such expressions, can be 
found in [3D] . For the sake of clarity, we report below the correspondence between the LCDAs used in this paper and 
those in [5U] : 



V2-K 



VP V^-K ' 

",^11 ^i-K ; 



*4:/f ; 



<t>3K -> ^''"^ ; 



A-K 
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